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ABSTRACT 

We report the detection of water maser emission from four nearby galaxies hosting ultradense H II (UDHII) 
regions, He 2-10, the Antennae galaxies (NGC 4038/4039), NGC 4214, and NGC 5253, with the Green Bank 
Telescope. Our detection rate is 100%, and all of these H2O "kilomasers" (Lh 2 o < 10 L ) are located toward 
regions of known star formation as traced by UDHII regions and bright 24 /im emission. Some of the newly 
discovered H2O masers have luminosities 1-2 orders of magnitude less than previous extragalactic studies and 
the same order of magnitude as those typical of Galactic massive star-forming regions. The unusual success of 
this minisurvey suggests that H 2 maser emission may be very common in starburst galaxies, and the paucity of 
detections to date is due to a lack of sufficient sensitivity. While the galaxy sample was selected by the presence 
of UDHII regions, and the UDHII regions lie within the telescope beam, in the absence of H2O spectral line 
maps the connection between H2O masers and UDHII regions has not yet been demonstrated. 

Subject headings: galaxies: interactions — galaxies: ISM — galaxies: star clusters — galaxies: starburst — 
masers — radio lines: galaxies 



1. INTRODUCTION 

Extreme star-forming environments known as "ultradense 
Hll regions" (UDHIIs) have been discovered in a num- 
ber of starburst ga l axies (e.g., Qj Cobulnickv & Johnson 19991: 
Turn er et al.1 120001: iTarchi et al.1 120001: Uohnson et al.1 120041: 
Reines et ani2008|) . UDHIIs are identified through their ther- 
mal mid-IR to radio emission (often being too deeply embed- 
ded in their birth material to observe at optical wavelengths). 
Similarly, dense radio Hll regions are known to exist on a 
much smaller scale around individual massive stars in the 
Galaxy (i.e., ultracompact Hll regions [UCHIIs]; Wood & 
Church well 1989). However, UDHIIs are vastly scaled up 
from individual UCHIIs, and they can contain the equivalent 
of a thousand or more embedded O-type stars. 

Given the tremendous number and density of young mas- 
sive stars in UDHIIs, it is logical to ask whether typi- 
cal signposts of Galactic massive star forma t ion ar e asso- 
ciated with these objects. IChurchwell et al.l (Il990h found 
that ~ 70% of Galactic UCHIIs are associated with H 2 
(22.235 GHz) masers. Subsequent interferometric follow- 
up showed that while some H2O masers are coincident with 
the UCHII regions themselves, many are instead associated 
with younger members of the forming m assive protocluster 
of which the UCHII is a m ember (e.g.. iTofani et al.l fl995b 
iHofner & Churchweil 1 19961) . These studies imply that the 
conditions required for H2O masers in > 10 7 cm -3 and T > 
400 K) are fairly persistent during the early stages of massive 
star cluster evolution. If extragalactic UDHIIs are composed 
of thousands of UCHIIs, then it is reasonable to expect H2O 
masers to be associated with a large fraction of the embedded 
massive stars in these giant clusters. 

The term "kilomasers" has been coined to describe extra- 
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galactic 22 GHz H2O masers with luminosities comparable to 
the brightest Galactic star formation-associated H2O masers 
(e.g., W49N; L H2 o ~ 1 L ). The luminosities of H 2 kilo- 
masers (Lh 2 o < 10 L ) are much lower than the more widely 
studied "megamasers," which are associated with the nuclear 
regions of AGN (i.e., NGC 4258) a nd have been observed 
with luminosities up to 10 4 L ( Barvain is & Antonuccil 
120051) . H2O kilomasers have been associated with both AGN 
and star formation activity; this luminosity regime includes 
the tails (high and low) of both maser populations. Although 
H2O kilomasers can be used to help pinpoint sites of active 
star formation, there have been few searches with the sen- 
sitivity necessary to detect these faint masers at extragalac- 
tic distances. To date six H2O kilomasers have been un- 
ambiguously associated with regions of known star forma- 
tion in the LMC, M82, IC 342 IC 10 M33 TIC 1331, an d 
NGC 2146 fWhiteoak & GardnedfT986l:lTarchi et al.l l2002albl: 
iHenkel et all 120051; ICastangia et al.l I2008L and references 
therein). The six remaining known H2O kilomasers either are 
associated with AGN activity or have ambiguous provenance. 
For example, the H2O masers in NGC 253 appear to be as- 
sociated with a nuclear outflow ( Ho fner et al.ll2006l) and not 
current star formation. 

We have carried out a Green Bank Telescope 5 (GBT) search 
for kilomasers toward four nearby starburst galaxies (3 < D < 
20 Mpc), known to host UDHIIs, down to a sensitivity level 
consistent with strong Galactic UCHII region H2O masers. 
Positive detections were found for all four galaxies as de- 
scribed in detail below. The unusual success of this minisur- 
vey is due to the fact that our observations are more than an 
order of magnitude more sensitive than the majority of previ- 
ous single dish surveys for either kilomasers or megamasers. 
This result suggests that H2O maser emission may be very 
common in starburst galaxies. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the 6i6-523 22.23508 GHz ortho-water maser 
line in He 2-10, the Antennae (NGC 4038/4039), NGC 4214, 

5 The National Radio Astronomy Observatory is a facility of the National 
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TABLE 1 
Journal of GBT Observations 



Galaxy 


Other Name 


Pointing 


Center 




D 


^int 


AVchan 


rms 


Observing 


UT Date 






a (J2000) 


5 (J2000) 


(Jy) 


(Mpc) 


(h) 


(km s" 1 ) 


(mJy) 


Mode 




He 2-10 




08 36 15.3 


-26 24 30 


26.4 


10.5(0.7) a 


2.2 


3.3 


1.2 


Nod 


2004 May 3, 4 


Antennae 


NGC 4038/4039 


12 01 55.0 


-18 52 58 


76.0 


20.0(1. 4) b 


2.1 


2.0 


0.9 


OffOn 


2004 May 4, 6 


NGC 4214 


UGC 7278 


12 15 39.2 


+36 19 40 


25.5 


2.94(0. 18) c 


1.1 


2.0 


1.3 


OnOff 


2004 May 3 


NGC 5253 


Haro 10 


13 39 55.8 


-31 38 20 


29.8 


3.33(0.17) d 


2.9 


0.66 


2.4 


Nod 


2004 May 3, 4 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. The 100 /im flux densities come from the IRAS 
Point Source Catalog. Ant indicates the total on-source integration time and AVdian is the velocity resolution of spectra plotted in Fig.0 

a For He 2-10 we adopt an Hi velocity-based Hubble flow distance that accounts for local peculiar velocities (Mould et al. 2000) and H Q - 72(3) km s -1 Mpc _1 . b The Antennae 
distance is controversial; we adopt the conventional value of 19.2 Mpc i Whitmore et al. 2005) updated to H Q - 72(3) km s -1 Mpc -1 . c The NGC 4214 distance is from the RG branch 



method i Maiz-Apellaniz et al. 2002). c 
range from 3.0 to 4.1 Mpc). 



The NGC 5253 distance is from Cepheids i Gibson et al. 2000), and the uncertainty is statistical only (other recent NGC 5253 Cepheid distances 



and NGC 5253 with the GBT on 2004 May 3-6 (Table 
Q]). Two 25 MHz offset IFs were observed simultaneously 
with 200 MHz bandpasses in two polarizations and 12.2 or 
24.4 kHz wide channels. A winking calibration diode and 
hourly atmospheric opacities were used for flux density cali- 
bration. Opacities ranged from 0.05 to 0.15 nepers but were 
typically below 0.10. The estimated uncertainty in the flux 
density calibration is 10%-15%. We used the 33" (FWHM) 
dual-beam K-band feed in a nodding mode for the compact 
galaxies (He 2-10 and NGC 5253) and a position-switched 
mode for the larger galaxies (the Antennae and NGC 4214). 
Bandpasses were Hanning smoothed to 48.8 kHz channels for 
a rest-frame velocity resolution of 0.66 km s -1 and Gaussian 
smoothed to the resolution listed in Table Q] to produce the 
final spectra shown in Figure [T] 

For all observations, records were individually calibrated 
and bandpasses flattened using the calibration diode and the 
corresponding off-source records, respectively. Scans and 
polarizations were subsequently averaged, and a fifth-order 
polynomial baseline was fitted and subtracted. The band- 
passes were rarely contaminated by radio frequency interfer- 
ence (RFI) but were carefully inspected and flagged. Base- 
lines were often not flat, but the atmospheric features in the 
bandpass were much broader than any possible maser lines. 
All spectra were carefully inspected by polarization, IF, and 
time subsets to confirm the validity of the detected lines. All 
data reduction was performed in GBTIDL, and all velocities 
are given in the heliocentric frame. 

3. RESULTS 

We have detected water maser emission in all four galax- 
ies observed. Measured line properties are listed in Table [2 
and further information for each galaxy is presented below. In 
all cases, given the complex velocity structure of the galaxies 
in this sample (likely due to the combined effect of large and 
small mergers and outflows), we might expect to find a corre- 
spondingly complex velocity structure connected to the maser 
emission. The precise association of the maser emission with 
specific regions and environments in these galaxies will only 
be possible with future high-resolut i on em ission maps. 

He 2-10: iKobulnickv & Johnson! (Il999h identified five dis- 
tinct UDHII regions within the blue compact dwarf (BCD) 
galaxy He 2-10, which is located at a distance of 10.5 Mpc. 
He 2-10 has roughly solar metallicity (12+log(0/H)=8.93; 
IVacca & Contilll992|) implying that molecules may be rela- 
tively abundant. The He 2-10 UDHII regions are associated 
with bright 24 /im emission, and the main body of the BCD 
fits completely within the 33" GBT beam as shown in Fig- 
ure [2]2, so at the current resolution it is not possible to pin- 
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Fig. 1. — Water maser spectra of He 2-10, the Antennae galaxies (NGC 
4038/4039), NGC 4214, and NGC 5253. Velocity resolutions (channel 
widths) are listed in Tabled Velocities are heliocentric. The bold horizontal 
bars indicate the velocity range of CO (see Tableland Numbers below 
the spectra label the lines enumerated in Tabled 

point where the H2O emission originates. The velocity extent 
of CO(l-O) emission in He 2-10 observed with OVRO is 800- 
940 km s~ l (IKobulnickv et al.ll 19951) . iMohan et al] (12001 ) find 
a H92a li ne velocity peak of ~ 950 km s -1 , with a broad blue 
wing, and lHenry et al.l (120071) find narrow Brackett line fea- 
tures in the range ~ 860-900 km s -1 . The maser velocity 
presented here of 887 km s -1 is roughly consistent with all of 
these studies. 

The Antennae: At an adopted distance of 20 Mpc, 6 
the interacting galaxy pair known as the Antennae 
(NGC 4038/NGC 4039) is the most distant source in our sur- 
vey. As shown in Figured, the single GBT pointing (33" 
beam) was centered on a 24 /im-bright region of copious 
star formation within the "interaction region" (IAR), offset 
from both interacting galaxies' nuclei. This region is also 
coincident with sever al CO-identified "s upergiant molecular 
complexes" (SGMCsi lWilson et alJ2000h a nd UDHII regions 
identified by VLA cm wavelength emission (Neff&Ulvestad 
2000). The broad component of maser emission that peaks at 
1515 km s -1 is consistent with velocities mea sured from ion- 
ized gas studies, such as Gilbert & Grahi3 d2007l) . as well 



6 Recent work by Saviane et al. ( 2004, 2008) suggests that the Antennae 
distance may be as close as 13.3 Mpc. 
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as that of CO toward the IAR (1350- 1700 km s" 1 ; see, e.g., 
IWilson et alJl2000t ISchulz et al.ll2QQ7h . However, the second 
maser component at 1670 km s -1 is on the high end of the CO 
gas velocities observed within the IAR, and may instead be as- 
sociated with a "bridge" of CO at 1650-1700 km s -1 that ap- 
pears to link the IAR t o the N GC 4039 nucleus in the channel 
maps of ISchulz etail (120071) . Emission from the NGC 4039 
nuclear region itself is unlikely because it lies near the first 
null of the GBT's beam; thus, Antennae feature 2 is unlikely 
to be associated with nuclear activity. 

NGC 4214: NGC 4214 is a BCD at a distance of 
2.9 Mpc with a subsolar m etallicity (12+log(0/H)=8.28; 
iKobulnickv & SkiUmanl[l996l) . As shown in Figure Eh H 2 
maser emission is detected toward a region of bright 24 /im 
emission near the center of the galaxy. The 24 /im and H2O 
emission are also coincident wi th several unreso lved ther- 
mal ra dio sources identifi ed by iBeck et aD (|2000|) . Using 
OVRO. IWalter et al.l (1200 lh found that the CO(1-0) emission 
in this galaxy is concentrated into three main regions, the cen- 
tralmost of which is coincident with the GBT pointing and 
the bright 24 /im emission. This CO emission has a peak 
V he i(CO)= 291.4 km s" 1 and a velocity extent of 287-308 
km s -1 , in good agreement with Vh e i(H20) = 301.6 km s _1 . 

NGC 5253: NGC 5253 is a dwarf starburst galaxy 
at a distance of 3.3 Mpc with subsolar metallicity 
(12+log(0/H)=8.23;|MMtin(i997|). Using OVRO. lMeier et al.l 
(2002) found that the CO(2-l) emission is concentrated to- 
ward the optically dark dust lane with a velocity range of 350- 
450 km s -1 . The GBT beam was cent ered on the UDHII re- 
gion discovered by Turner et al. ( 2000), and a region of bright 
24 /im emission (see Fig.[2ji). This region is also coincident 
with the weak CO(2-l) molecular cloud "D" with Vh e 1 (CO)= 
394 km s" 1 and Vfwhm=15 km s" 1 (iMeier et al.ll2002l) . How- 
ever, the GBT beam also encompasses parts of CO cloud "C" 
(Vhei=418 km s -1 ) to the southeast and cloud "A" (Vh e i=363 
km s -1 ) to the east. Thus, the H2O maser emission at 366 and 
375 km s -1 is significantly blueward of the CO gas toward 
the UDHII region and most closely agrees with the velocity 
of clo ud "A." Using the radio recombinations lines H53a and 
U92a iMohan et al.l (1200 lh and iRodnguez-Rico et al.l (l2007h 
measure peak velocities of ~ 400 km s" 1 , both with broad 
wings. These velocities are somewhat higher than, but con- 
sistent with, the maser velocities. Given that the ionized 
gas in N GC 5253 is known to hav e a sinusoidal velocity 
structure (|L6pez- Sanchez et all F2007) with an amplitude of 
~ 20 km s" 1 , it is not surprising that we see maser lines with 
velocity separations within this range. 

4. DISCUSSION 

Could the strong Galactic H2O maser source W49N be 
detected in our sample galaxies? Using a fiducial dis- 
tance of 11.4 kpc and a 10 kJy peak flux density (see e.g., 
iLiljestrom et al.lll989|) . W49N is marginally detectable in the 
Antennae (3 mJy) but easily detectable in the other sample 
galaxies (12-150 mJy). By luminosity (~ 1 L ), W49N 
should be detectable in NGC 4214 and NGC 5253 and possi- 
bly in He 2-10, but not in the Antennae. It is remarkable that 
we can now detect sub-W49N H2O masers in external galax- 
ies, which means that extragalactic H2O kilomaser observa- 
tions may now be interpretable in a Galactic star formation 
context. 

Two of our new H2O maser detections are extremely low 
luminosity; there is only one other system known with a total 
Lh 2 o < 0.1 L , and this was a brief flaring episode in IC 342 



(iTarchi et al.l l2002al: iHenkel et all 120051) . The detection of 
deci-W49N emission in a larger sample of galaxies similar to 
NGC 4214 and NGC 5253 will provide a measurement of the 
H2O luminosity functio n at least an order o f magnitude below 
current limits on L Hl o (IHenkel et al.l l2005) and will genuinely 
explore the link between H2O masers and star formation in 
galaxies as a function of environment, metallicity, and age. 

Our minisurvey sample has a 100% detection rate, and ex- 
pands by 2/3 the sample of known H2O kilomaser galaxies 
associated with star formation. The GBT beam excludes the 
nuclei in the Antennae, so the association of H2O maser emis- 
sion with star formation is likely in this galaxy merger. More- 
over, the three dwarf galaxies in our sample show no evidence 
of AGN activity. Both He 2-10 and NGC 5253 have been 
observed with the very long base line interferometry (V LB I) 
High-Sensitivity Array (HSA) bv lUlvestad et al.l (120071) . and 
neither has pointlike radio sources down to a detection level 
of la ~ 21 /iJy, consistent with no AGN acti vity. The mid- 
infrar ed spectrum of NGC 4214 was studied by Sat yapal et al.l 
(|2008|) . and they find no evidence of high-ionization lines that 
would be suggestive of an AGN. Indeed, given the associa- 
tion b etween the presence of AGN and host galaxy bulge mass 
(e.g., KauffmaZetiDl2003), it would be surprising if any of 
these three dwarf galaxies harbored an AGN. Without H2O 
spectral line maps, however, the connection between H2O and 
UDHII regions has not yet been demonstrated. 

All four detected masers lie > 1 dex below the log(L H2 o) = 
log(L FrR )-9.5 relat ionship obtained for Galactic H2O masers 
by felli et al. (Il992h . That is, these masers are "sublumi- 
nous" in the H2O line given their host galaxies' far-IR lu- 
minosities. They are, however, in good agreement with the 
\ogL H2 o = logLf/R - 10 relation ship for H2O kilomasers ob- 
tained by lCastangia et al.l (|2008|) . The origin of this difference 
between kilomasers and Galactic masers is unknown, but is 
likely due to spatially unresolved far-IR (IRAS) observations 
of galaxies that include IR emission unrelated to the sites of 
maser emission (Galactic sites of star formation were much 
better resolved). The difference between the H2O-FIR rela- 
tion of kilomasers and Galactic masers may thus not be phys- 
ical at all (unlike H2O megamasers, which represent a sub- 
stantially different physical process from Galactic and kilo- 
masers). 

5. CONCLUSIONS 

We have identified H2O masers in four galaxies hosting ul- 
tradense Hll regions. The high detection rate and the low 
luminosity of the H2O lines in our minisurvey indicate that 
previous water maser surveys of nearby IR-bright galaxies 
lacked the sensitivity to detect most H2O kilomasers associ- 
ated with bursts of star formation. We suggest that sensitive 
future surveys will produce a wealth of new H2O maser de- 
tections, enabling detailed studies of intense star formation in 
other galaxies. Future interferometric imaging of the maser 
lines detected with the GBT are needed to conclusively deter- 
mine whether the the maser emission is connected to massive 
star formation in these galaxies. Indeed, although each of the 
GBT-detected masers are likely composed of numerous indi- 
vidual maser spots, if any of them remains relatively bright at 
very high spatial resolution (rather than segregating into many 
weak spots), it may be possible to carry out VLBI astrometry 
in the future as has been recently done for several Lo cal Group 
galaxies hosting H2O masers dBrunthaler et al.ll2007|) . 

This research has made use of the NASA/IPAC Extragalac- 
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Fig. 2.— Three-color archival Spitzer images of (a) He 2-10, (b) the Antennae, (c) NGC 4214, and (d) NGC 5253 with RGB mapped to 24, 8, and 3.6 /mi. In 
these images red colors pinpoint the most active regions of star formation. The locations of the 33" GBT beam are indicated with green circles; the positions of 
the two Antennae nuclei are also indicated. The GBT beam corresponds to a physical scale in these galaxies of 1.7, 3.2, 0.47, and 0.53 kpc, respectively. 



TABLE 2 
Water Maser Line Properties 



Galaxy 


Line 


Vh 2 o 


Vco 


Speak 


AVpWHM 


JSdv 


Lh 2 o 






(km s" 1 ) 


(km s" 1 ) 


(mJy) 


(km s" 1 ) 


(mJy km s 1 ) 


(Lq) 


He 2-10 


1 


887(6) 


800-940 


2.4(0.3) 


115(14) 


267(10) 


0.68(9) 


Antennae 


1 


1514.7(4) 




12.4(0.9) 


24(2) 


704(10) 






2 


1670.4(2) 




14.5(0.9) 


7(1) 


172(4) 






Sum 




1350-1700 






892(12) 


8.2(1.1) 


NGC 4214 


1 


301.6(8) 


287-308 


6.7(0.6) 


21(2) 


141(6) 


0.028(4) 


NGC 5253 


1 


366.0(1) 




18.1(2.6) 


1.2(0.2) 


29(8) 






2 


375.1(2) 




13.2(1.6) 


3.3(0.5) 


47(9) 






Sum 




350-450 






83(13) 


0.021(4) 



Note. — The line properties were determined from Gaussian fits, with the exception of the Antennae, which were obtained directly from the data because the line profiles are non- 
Gaussian. The peak channel in the He 2-10 spectrum is 4.5 a alone; the gaussian-fit peak is ~8 a. Isotropic line luminosities were computed from Lh 2 o = (0.023 L©) x D 2 x J Sdv, 
where D is in Mpc and J Sdv is in Jy km s -1 I Henkel et al. 2005). Vco is the range of CO line velocities previously observed within the GBT beam; references are provided in § 3. 
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